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It has been recently suggested that the magneto-rotational instability (MRI) is a limiting 
case of the magneto-elliptic instability (MEI) . This limit is obtained for horizontal modes 
in the presence of rotation and an external vertical magnetic field, when the aspect ratio 
of the elliptic streamlines tends to infinite. In this paper we unveil the link between these 
previously unconnected mechanisms, explaining both the MEI and the MRI as different 
(-H ■ manifestations of the same Magneto- Elliptic- Rotational Instability (MERI). The growth 

I rates are found and the influence of the magnetic and rotational effects is explained, 

in particular the effect of the magnetic field on the range of negative Rossby numbers 
at which the horizontal instability is excited. Furthermore, we show how the horizontal 
rotational MEI in the rotating shear flow limit links to the MRI by the use of the local 
shearing box model, typically used in the study of accretion discs. In such limit the growth 
rates of the two instability types coincide for any power-type background angular velocity 
radial profile with negative exponent corresponding to the value of the Rossby number 
of the rotating shear flow. The MRI requirement for instability is that the background 
angular velocity profile is a decreasing function of the distance from the centre of the 
disk which corresponds to the horizontal rotational MEI requirement of negative Rossby 
numbers. Finally a physical interpretation of the horizontal instability, based on a balance 
between the strain, the Lorentz force and the Coriolis force is given. 
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1. Introduction 

The stability of vortices is a problem of paramount importance in fluid mechanics. 
Considering that turbulence consists of vorticity blobs, vortices are the fundamental unit 
of turbulent flow. Unveiling the mechanism that renders them unstable should therefore 
provide vital insights into the nature of turbulence itself. Phenomenologically, turbulence 
can be described as a series of bifurcations, starting with a primary instability that con- 
verts shear into vorticity, creating vortices. This is followed by another bifurcation, a 
secondary instability, to break these vortices into lesser vortical structures. These in turn 
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shall experience a sequence of "inertial instabilities", leading to a cascade. Though the 
Kelvin-Helmholtz instability and the Ray leigh- Taylor instability are well established as 
examples of primary instabilities, the highly successful theory of the turbulent cascade 
put forth by Kolmogorov (1941) rested on a heuristic picture of secondary instability, 
established by early experiments (e.g., Taylor 1923). It was not until the works of Picr- 
rehumbert (1986) and Bayly (1986) that the elliptic instability was introduced as a 
mechanism for the secondary instability. A fluid in rigid rotation supports a spectrum of 
stable inertial waves, the simplest case being circularly polarized transverse plane waves 
oscillating at twice the frequency of the base flow. Strain is introduced when the stream- 
lines pass from circular to elliptical, and some modes find resonance with the strain field, 
leading to de-stabilization. The nonlinear evolution and saturation of the elliptic insta- 
bility was studied numerically by Schaeffer & Le Dizes (2010) and experimentally by 
Eloy et al. (2003). Herreman et al. (2010) generalized the linear and weakly nonlinear 
theories for the elliptic instability (WalefFc 1989) to the magnctohydrodynamics case and 
conducted experiments to explain some aspects of the nonlinear transition process. 

The seminal work of Picrrehumbert (1986) and Bayly (1986) was followed by studies 
that considered the incorporation of more physics into the problem. Miyakazi (1993) 
included background rotation, unveiling a strong two-dimensional mode, for which the 
wavevector and the vortex axis are parallel. For incompressible flow such modes are 
solenoidal, so the motion itself is horizontal to the vortex axis. This horizontal instability 
is thus the two-dimensional mode of what could be called an elliptic-rotational instability 
(ERI). This instability is not of resonant nature, but centrifugal, appearing as exponential 
growth of epicyclic disturbances. This behaviour invites a connection with the Raylcigh 
(centrifugal) instability, and indeed the mechanism is similar (Lesur & Papaloiozou 2009), 
suggesting that the Rayleigh instability is a limit of the ERI. 

If the elliptic instability is such a fundamental process for understanding hydrodynami- 
cal turbulence, one can expect that the magneto-elliptic instability (MEI) should likewise 
provide a similar framework when it comes to MHD turbulence. In the geophysical con- 
text, Kerswell (1993) studied the influence of a toroidal magnetic field and stratification, 
both axial and radial, on the stability of elliptic flow in uniform rotation. In his con- 
figuration, axial stratification and the magnetic field were invariably stabilizing, though 
radial stratification enhanced the elliptical instability. Lebovitz & Zweibel (2004) stud- 
ied the problem of stability of axially magnetized elliptic streamlines, finding that the 
MEI has new families of destabilizing resonances, corresponding to Alfven waves. Just as 
Miyakazi (1993) extended the work of Bayly (1986) to include background rotation, the 
work of Lebovitz & Zwebeil (2004) was similarly extended by Mizerski & Bajer (2009), 
who analysed the joint effect of the magnetic and rotational effects on the stability of the 
Eulcr flow with elliptical streamlines, by the use of numerical and asymptotic analytical 
methods in the limit of small ellipticity of the flow. Their comprehensive study, based on 
the analytic technique developed by Lebovitz & Zweibel (2004) included both the modes 
excited via a resonance mechanism and the horizontal modes, that propagate along the 
magnetic field lines. The case of the horizontal instability was then further investigated 
in the non-magnetic case by Lesur & Papaloizou (2009) and in the magnetized case by 
Bajer & Mizerski (2011). However, the problem of the MEI in the presence of rotation 
is conceptually difficult and there are aspects which still require clarification and further 
insight is necessary. 

One such insight may come from astrophysics, which also provides the motivation for 
the present study. When modeling the problem of magnetized vortices in protoplane- 
tary disks, Lyra & Klahr (2011) found that the horizontal magneto- elliptic mode, when 
taken to the limit of infinite ellipticity, yields the same growth rates and range of unsta- 
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ble wavenumbers as the well known magneto-rotational instability (MRI, Velikliov 1959, 
Chandrasekhar 1960). The MRI is an instability of magnetized Eulcr flow in circular 
differential rotation, which renders the motion unstable when the angular velocity de- 
creases outwards. Since Keplerian flow is an Euler flow in circular differential rotation, 
and magnetic fields are ubiquitous in astrophysics, the MRI is the best candidate to gen- 
erate turbulence in accretion disks (Balbus & Hawley 1991, 1998, Armitage 1998, Hawley 
2000, Fromang & Nelson 2006, Lyra et al. 2008, Flock et al. 2011). 

The result of Lyra & Klahr (2011) tantalizingly suggests that the MRI is a limiting case 
of the MEL In this paper we explore this interesting prospect in more detail, unveiling the 
connection between these two previously unconnected instabilities. In a way, the present 
study is an attempt at unification, explaining both instabilities as different manifestations 
of the same Magneto- Elliptic- Rotational Instability (MERI). 

The paper is structured as follows. In the next section we briefly review the horizontal 
rotational MEI presenting a somewhat different approach from Bajer & Mizerski (2011) 
and show that the external magnetic field has a profound effect on the range of values 
of the Rossby number for which the horizontal rotational MEI is excited. In Sect. 3 we 
further develop the idea presented in Lyra & Klahr (2011) and demonstrate, by the use 
of the so-called local shearing box model (cf. Hawley et al. 1995) how the horizontal 
rotational MEI, in the limit of rotating shear flow can be related to the MRI. Next we 
present a physical interpretation of the horizontal instability based on a balance between 
strain generated by the basic flow and the Coriolis and Lorentz forces. Conclusions are 
presented in Sect. 5. 



2. The calculation of growth rates of the horizontal rotational MEI 

Let us consider the Euler flow with elliptic streamlines, 

uo=j[-{l + e)yAl-e)x] (2.1) 

rotating with the angular velocity n = ile^, where 276^ is its uniform vorticity and 
< e < 1 is the strain. We now briefly review the stability properties of the flow (|2.H) in 
the presence of uniform vertical field Bo = Boe^ with respect to the horizontal modes, 
which are defined by the fact that they propagate only in the vertical direction, i.e. the 
wave vector of the perturbation has only the 'z '-component. The solenoidal constraints for 
the velocity and magnetic fields then imply that the perturbations have only horizontal 
components. Thus the velocity and magnetic field perturbations have the form 

v = ve"=^ and b = be"=^ (2.2) 

respectively. The equations to solve are the Euler and induction equations 

dm + vjdv, = --d,P - 20^jv.j + —BkdkB^, (2.3) 

dtB,+v,djB, = BjdjB,, (2.4) 

where the magnetic pressure and centrifugal potential were incorporated in the pressure 
P and Oik = £ijk^j- We apply the scaling hi — bi/y/^oP so that, in principle, we are 
solving for the Alfven velocity and after linearisation the evolution equations for the 







K. A. Mizerski and W. Lyra 






Horizontal MERI - Unstable range 




2.0 


/7=0 /' ' 

■ h=2 i ' 

- h=4 i I 


1 ' 




1 .5 


i 1 
/ 1 
' 1 
/ 1 

/ 1 
/ ' 


Ro"' - -lr/4 . 




1.0 












'o 
cc 


















0.5 




instability 




0.0 




Rq-' ^ -(h- + I)/2 _ 






■ . , , 1 . . , 1 . . . 1 , , . 1 . , , 





-10 -8 -6 -4 -2 

Ro-' 



Figure 1. The square of the function ^+ as in Eq. (|2.7[l . that defines the unstable horizontal 
modes of the magneto-elliptic-rotational instability. Instability exists when \ip+\ < e, where 
^ e ^ 1 is a measure of eccentricity. The instability exists only for negative Rossby numbers, 
i.e., anticyclonic rotation. The unstable range of wavenumbers is Ro~^ < —h? /A and unbounded. 
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(2.5) 



where t = t/j and the upper dot denotes a derivative with respect to r, Ro = "f is the 
Rossby number and h = kBo/jy'iJ-op. The eigenvalues of the matrix in (|2.5p which we 
denote by Cj, j = 1, 2, 3, 4 are the complex growth rates of the perturbations. They are 



(y-i h ^ -'^IPl, (74/7 = -y^^^ -V': 

where 

V'± = Ro"^ ± \/(Ro"^ + i)%;i2^ 

Note that for all /i > and Ro e M 



(2.6a) 
(2.66) 

(2.7) 

(2.8) 



which means that 0-3 and 0-4 are imaginary, thus produce only oscillations. The basic 
state is unstable with respect to perturbations propagating along the 'z' axis if and only 
if there exist Ro and h such that 

< 1 and e^lV'+l- (2.9) 
The function ?/;_|_ (Ro~^; /i), plotted in Fig. [TJ possesses the following properties 



/,2 

V'+ ^ 1 ; for Ro-i ^ - — 
t/j-h = 1 ; for Ro"^ = -h^jA 



(2.10a) 
(2.106) 
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(2.10c) 
(2.10rf) 

(2.11) 



lim Tp^ (Ro ^;fij = 

Ro~ ^ ^ — oo 

lim Ro"^ {l - [l + RO + C' (Ro^)]} = -1. 

Ro ^ — )■ — oo 

Therefore the instability is only possible if Ro^^ < —h? /A since only then, for a certain 
ellipticity e ^ jV'+l the growth rate ai can be real and positive. Moreover, equation 
(|2.12p . in particular, means that as Ro~^ tends to — oo, for some Ro~^ < ~{h? + l)/2 
the function ?/;+(R-0~"'^; h) must have an inflection point. The limit /i — > is singular in 
the sense that the two curves '0^(Ro~^; h) and '!/'?.(R.o~^; h) merge into a parabola (and 
a straight line parallel to the Ro~^ axis of constant value 1), touching each other at the 



inflection point of the function -0^ (Ro ; h) at Ro 



-1, = "0- = 1- This means 



that for /i = the horizontal instability is possible only for a bounded set of values 
Ro"""^ G (—1,0). The situation significantly changes immediately after switching on even 
a very weak magnetic field since then the domain in which the horizontal instability is 
possible becomes unbounded, i.e. Ro~^ S (— oo, —h? /A). As the instability exists only for 
negative Rossby numbers, only anticyclonic vortices are unstable. 

As Lyra & Klahr (2011) have pointed out, there is instability in the limit of a vanishing 
vortex in a shearing box (cf. Hawley et al. 1995). In that case, the Rossby number remains 
finite owing to the vorticity of the Keplerian shear, and the case is equivalent to a Kida 
vortex (Kida 1981) of infinite aspect ratio 



1 



and Ro = — 
4 



(2.12) 



In the limit e = \ only the horizontal instability persists (the purely horizontal modes 
are the most unstable ones and the modes destabilised via resonances are suppressed; cf. 
Mizerski & Bajer 2009, 2011), for which the dispersion relation is (cf. equation (|2.6b ')') 



= £ — 



r 



Ro~ 




(2.13) 



where va = Bo/y/JIop is the Alfven speed. For Ro < 0, also substituting 7 = Roil and 
q = kvA/^, we obtain 

n 2 



= - Ro" 

7 



(Ro + 1) 



g2 - 1 



and applying the limits defined by (|2.12p 

72 ^ 9 



1-2^^-1 



(2.14) 



(2.15) 



Thus the growth rate is real and positive for < g < \/3, and the maximal value 
cr = I7I = (3/4) \Q\ is achieved for q = \/l5/4 ~ 0.9682. This is precisely the same 
dispersion relation obtained for the MRI by Balbus & Hawley (1991, see also Hawley & 
Balbus 1991). 
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3. Why do we see the MRI in a rotating shear flow? 

The MRI was originally formulated for axisymmetric flows with circular streamlines 
and the condition for MRI to operate is that the angular velocity decreases with radius. 
However, as it turns out, the MRI also seems to operate in non-axisymmetric flows and 
appears as a limiting case of the MEI in the presence of rotation, in the limit when 
the latter becomes a rotating shear flow (Lyra & Klahr 2011). Let us now assume that 
the Euler flow with elliptic streamlines given in (|2.ip . rotates with the angular velocity 
ft = rte^ about an axis parallel the 'z' axis and located at distance R from the centre of 
the elliptical vortex. Let us also consider the inertial (non-rotating) system of reference 
{x',y',z'), with the origin x' = 0, y' = at the axis of rotation (see Fig. In this frame 
of reference the vector of position of the centre of the elliptical vortex is 



R' = R[-sinnt, cos fit] , 
where the prime indicates representation in the non-rotating frame and 

dR' 

"dT 



= rj X R', 



(3.1) 



(3.2) 



where the time derivative is taken within the non-rotating frame, i.e. holding the unit 
vectors e' and e' constant. It was chosen that 



Using the transformations 

cos fit sin fit 
— sin fit cos fit 

and 



X 









R' (< = 


= 0) = 




x' 









v' 




R 



cos flte'^ - 
— sin flte' 



sin flte'y 
f cos flte[, 



flxr' 
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rjxR' = u + r2xr' 



(3.3) 



(3.4) 



(3.5) 



where r' = [x' ,y' , z']^ , we obtain the following expression for the basic elliptic velocity 
field in the inertial frame 



'-'qx = -7 W + £ {y'cos2flt - x' sin 2m)] + (l + e) R cos fit - fly' , 



(3.6) 



u'^y =,^[x' ~e {x' cos 2flt + y' sin 2^)] + {I + e) R sin fit + fix' . (3.7) 

Using the cylindrical coordinates {r' ,ip' ,z) we may write the angular velocity of the flow 
in the inertial frame as 



r' 



u'^yx' u'oxU' 

r,Jl 



R 



= f7 <^ 1 + Ro [1 - £ cos (2(p' - 2flt)\ - Ro (1 + e) — sin (^' - fit) 



(3.8) 

which, clearly, is a function of r' (as long as i? 7^ 0) and ^p' (in the above Ro = 7/^7 is the 
Rossby number). Lyra & Klahr (2011) have shown numerically that in the rotating shear 
flow limit defined in (|2.12p the rotational MEI growth rate (i.e. the horizontal instability 
growth rate) profile with the parameter q indeed matches the MRI profile. In such limit 
the expression p.Sp yields 



(3.9) 



!f2^ ^ J7 (1 + ^ cos (2(^' - 2flt) + -- sin U' - flt)\ 
r' 1^ 4 4 2 r' J 



therefore the angular velocity in the inertial frame, in the rotating shear flow limit, is 
depended on the distance from the origin as long as i? 7^ 0. Moreover, the flow in the 




X 

Figure 2. A schematic picture of the rotating elhptic flow. 



inertial frame is not axially symmetric, but nevertheless the resemblance to the MRI 
is observed. Furthermore the linear stability analysis which leads to the expressions for 
the growth rate (|2.13p and thus (|2.14|) is independent of the value of R which means 
that even for i? = 0, i.e. when the angular velocity is independent of r' we obtain the 
MRI-like results in the rotating shear flow limit of the rotational MEL Since the Rossby 
number must be negative for the instability to settle in, it seems that it is crucial that 
the MRI requirement for the angular velocity to be a decreasing function of r' must 
be met in the neighborhood of the centre of the rotating elliptic vortex, i.e. for x = 0, 
2/ = 0. In principle for a given value of the Rossby number the angular velocity (|3.9p is 
a decreasing function of r' only in a half-plane and an increasing function of r' in the 
other half-plane, depending on the signs of Ro and sin ((^s' — Vtt). However, as said, the 
instability occurs only for Ro < 0, in which case the angular velocity decreases with r' 
for Vlt<if'<Ti + fit, and the steepest decrease is observed a.t ip' = tt/2 + fit, which is 
the azimuthal coordinate of the centre of the elliptic vortex. 

The explanation is based on the local shearing sheet model for accretion disks, since in 
this model the shearing flow is simply a local representation of the background rotation 
of the axisymmetric vortex with angular velocity decreasing with r' . In other words for 
an axisymmetric vortex with the angular velocity depended on r' as ^{r') ^ l/r'*' and 
thus gravity as g ^ r'", with n = —2s + 1, in the vicinity of r' = i? we get 

^(r) ^ 1 ^-—s~—-sr^—-, (3.10) 

and from p.8p . in the vicinity of r' = R', i.e. for r' = R + r and Lp' = Tr/2 + fit + Lp, we 
obtain 

^ = 17(l + Ro + Rocos(2v3)-2Rocos(/3 + 2Ro— cos(p|. (3.11) 
r' L R J 

According to the shearing sheet approximation we now substitute in (|3.1ip for 51 = 1 /i?* 
and since cosip 1 + 0(^9^) the expression p. lip up to order one in the distance from 
the point r' = R' (thus in the region of validity of the shearing sheet approximation) 
takes the form 



which agrees with p.lOp for 



Ro=-^ = ^. (3.13) 
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The Kida (1981) solution for a stationary elliptic vortex embedded in a shear flow yields 



which in the shear flow limit e = 1 agrees with p.l3p . The horizontal rotational MEI in 
the presence of a magnetic field operates for all negative Rossby numbers and according 
to (|3.13p this corresponds to s > 0, i.e. the angular velocity ri(r') decreasing with r', 
which is the requirement for the MRI. It is now clear, that since the growth rate of the 
horizontal instability in a rotating shear flow (|2.13p docs not depend on position, and at 
the centre of the elliptical vortex r' — R' the rotating shear flow models locally the flow 
of an axisymmctric disk with angular velocity decreasing with r', the growth rate must 
correspond to the growth rate of the MRI taken at r' = R', since in the short-wavelength 
limit analysed by Balbus & Hawley (1991) the MRI criterion is local. Nevertheless it 
remains an interesting feature of the rotating shear flow limit that despite the fact that 
the flow is not axially symmetric and for i? = can even be independent of r' the stability 
characteristics are the same as of the classic axially symmetric, short-wavelength MRI. 

It is also important to realise that, in general, the horizontal instability persists even 
without the presence of the magnetic field. The MRI operates only if the magnetic field 
is present. Furthermore, as already mentioned, even the r'-dependence of the angular 
velocity is not necessary for the development of the horizontal instability, since if i? = 
the expression (|3.8p becomes independent of r' and the growth rate (|2.13|) is not affected. 
Thus, in principle, it is the lack of axial symmetry that is crucial for the development 
of the horizontal rotational MEI. Nevertheless there exists a common physical ground 
for both the instability types, i.e. horizontal rotational MEI and MRI, since the rotating 
shear fiow limit of the rotational MEI, even though non-axially symmetric, exhibits the 
properties of the MRI. The common feature of both instability types is the shear which 
is present in the elliptic fiows as well as in axisymmctric flows with angular velocity 
dependent on radius and this feature accompanied by the rotational effects (local in MRI 
and global in rotational MEI) is behind the physical interpretation of the common limit 
for both instabilities discussed in here. This point is elaborated in the next Section. 

We have, therefore, demonstrated the direct link between the megncto-elliptic and 
magneto-rotational instabilities and in the non-magnetic case between the horizontal 
elliptic and centrifugal instabilities and the physical conditions under which previously 
unrelated instability types can be linked were explained in detail. The dispersion relations 
for horizontal MEI and MRI, likewise for horizontal EI and the centrifugal instability are 
the same in the limit e = 1, which suggests, that the physics of the horizontal magneto- 
elliptic instability is similar to the physics of the MRI. 

Let us make a final comment on the non-magnetic case. As mentioned the horizontal 
instability is present even if the magnetic field is not. The axisymmetric flows in the 
non-magnetic case are only subject to the centrifugal instability for which the Rayleigh's 
criterion says that instability can occur if the angular momentum is a decreasing function 
of r' in some region. According to p. 101) , p.l2[) and p.l3|) the angular momentum of an 
axisymmetric vortex decreases with r' if the Rossby number for the rotating shear flow 
satisfies Ro < — 1. This is an exact instability condition for the horizontal modes in the 
non-magnetic case (cf. Mizerski & Bajer 2009) for a rotating shear flow, i.e. for e = 1, 
since in such case the growth rate of the horizontal modes is (cf. (j2.13|) and p.l3|) ) 



s 



2 



(3.14) 



Ro 




CT^ = 2 (s - 2) . 



(3.15) 



On the other hand, the classic Rayleigh's criterion, formulated for systems with imper- 
meable boundaries at r' = Ri and r' = R2, for horizontal modes and ri(?'') ~ l/r'" 
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yields 



2R 



2-2s 



2-2s 



(3.16) 



s~l Rl-Rl ■ 

If the Rayleigh's criterion is made local, by assuming a narrow gap limit, i.e. R2 = Ri+6r 
and Sr/Ri <C 1 the above expression (|3.16p reduces to 

'^-2(s-2)n^ (i?i). 



(3.17) 



which agrees with (|3.15p . 

It seems in place to briefly comment on the case when the applied field in a disc is 
azimuthal, since the numerical and experimental works of HoUcrbach & Riidiger (2005) 
and Stefani et al. (2006) suggest a great importance of the azimuthal field for the re- 
duction of the instability thereshold. For the rotating shear flow Uq = —2jyex modeling 
locally the flow in a disc one could take Bq = (Bo + Joy)^x- Such stability problem with 
perturbations in the form of horizontal modes takes the form 



with by ~ and F 
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(3.18) 



Solving for the eigenvalue one obtaines 



— = -mo- 

T 



[1 + Ro- 



F^ 



(3.19) 



which leads to a conclusion that the instability is possible if either Ro < —2/(1 + VT+T^) 
or Ro > 2/(\/l + F^ — 1) is satisfied. Thus a new unstable branch appears for large and 
positive Rossby numbers, which has the property, that for weak field gradients, F ^ 1 
only strongly shearing discs are affected, Ro = —s/2 ^ 1. In principle the correspondence 
between the stability characteristics of the rotating shear flow with horizontal field and 
the local stability characteristics of a disc with azimuthal field should exist, however, the 
horizontal modes anaysed here are independent of y, thus axisymmetric in the global 
sense and as argued by Terquem & Papaloizou (1996) the most unstable modes in the 
case of azimuthal fields are the non-axisymmetric ones. Therefore the dynamics of discs 
with azimuthal basic fields should be dominated by the growth of non-axisymmetric 
modes not analysed in here. 



4. Physical interpretation of the horizontal rotational MEI 

The essence of the rotational MEI is the departure from axial symmetry, which im- 
poses strain on the perturbations. The typical model based on consideration of radial 
perturbations which preserve angular momentum, useful for physical interpretation of 
the MRI and the centrifugal instability is no longer applicable when the vortex is ellip- 
tical. A similarly simple physical model of the rotational MEI would have to take into 
account not only the elliptical shape of the vortex but also the action of the Coriolis 
force, which is essential for the instability mechanism. We present a somewhat simpler 
approach and provide the physical explanation of the horizontal rotational MEI based 
on the balance of forces acting on a general horizontal perturbation, highlighting the 
importance of strain. Let us write down the equations governing the evolution of the 
horizontal modes, for which the wave vector is vertical, k = fce^, and thus, because of 
the divergence free constraint, the perturbations posses only the horizontal components. 
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Figure 3. A sketch of the forces acting on fluid element in the presence of horizontal 

perturbation. 



The Navicr-Stokcs equation in the rotating frame and in the non-magnetic case is 



9v 
'dt 



+ (v • V) uo = -2J1 X V, 



(4.1) 



where Uq is the basic elliptic flow defined in (j2.1[) and v is the horizontal perturbation 
in the rotating frame. Note that the term describing advection of the perturbation by 
the elliptic flow, (uq ■ V)v and the pressure terms vanish in the case of purely horizontal 
perturbations and the only terms that remain are the terms describing the stretching 



of the perturbation by the basic flow and the Coriolis force. Introducing v 
T = jt we can rewrite the equation (|4.1|) in the following form, 



and 



_d_ 
d7 



= 2Ro" 



1) 



(4.2) 



where Ro = 7/17 is the Rossby number. The first term on the right hand side of (|4.2p 
corresponds only to local rotation and the second one to the local rotation and strain (in 
the following we will refer to the first term as the 'pure rotation term' and to the second 
one as the 'straining term'). We can easily solve the above equations (|4.2p to obtain 



1 



crr/7 



1 

+ 2 



- Vxo) < 



r/7 



(4.3) 



h = \ i^yo + XVm) e""/'' - ^ {-Vya + xv^o) 6"""/^ (4.4) 
where we have assumed that > (l + 2Ro^^)^ (thus the system is unstable) and 



X 



|'l + 2Ro"^) 



l + 2Ro" 



a 

1 



fl + 2Ro" 



(4.5) 



(see p.lSp 'l. The constants Vxo and VyQ are the initial values of the components of the 
perturbation. It can easily be seen that in the case when no strain is present, i.e. e = 0, 
the only term that remains in the equation (|4.2|) corresponds to local rotation and there 
is no possibility for instability (cf. the growth rate in (|4.5p ). If, on the other hand, the 
pure rotation term vanishes, i.e. Ro = —2 and only the straining term remains, the 
destabilisation is maximal. The unstable mode, which dominates the dynamics at large 
T has a direction defined by 



^ (unstable) 

^jy 

^(unstable) 



X- 



(4.6) 
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In general the straining term possesses not only the component parallel to the perturba- 
tion but also the perpendicular one (see figure 2) and they are given by 



Fl = (4.7) 



II _ 2eVxVy 

V 

vl — V. 



F^=e^ ^, 4.8 

V 

where v = -^^x + the magnitude of the pure rotation term, which is always 

perpendicular to the velocity vector of the perturbation is 

F2 = (l + 2Ro-i)'w2. (4.9) 

The parallel component fJ can either amplify the perturbation or damp it depending on 
the sign of VxVy, hence if the two perpendicular forces and Fn can not be in balance 
for any and Uy, the perturbation will not grow for all times, but it will oscillate. The 
perpendicular forces can only be in balance if 



2Ro"^) 



2 

> 1, (4.10) 



and then the perturbation grows unboundedly, amplified by the parallel component Fg 
(cf. Fig. [3]). The equation (|4.10p is solved by any perturbation satisfying (j4.6|) . which 
means that if the balance between the perpendicular forces can occur for some {j^""^''''^''^) 
and 1)^""*'*'''^''') ^ QT^y random perturbation will tend to the unstable mode defined by (j4.6p 
by rotating the perturbation vector towards the unstable direction (and only one di- 
rection is stable, for which the perpendicular forces are also in balance but the parallel 
force acts in opposite direction to the direction of the perturbation velocity vector damp- 
ing the perturbation, i.e. for ^'^^^^^^^ jy'^'^^'^^^) _ Furthermore, the equation (|4.10p 
means that destabilisation can only occur if > (l -t- 2Ro^^) , which is the instability 
requirement obtained in (|4.5p . 

In the rotating shear flow limit, e = 1, which as shown earlier corresponds locally to 
the centrifugal instability with horizontal modes, the unstable direction is determined 
by X = ■\/— 1/ (1 + Ro). Using the transformations from the previous section (j3.4p we 
can obtain the form of the unstable mode in the non-rotating frame, in the vicinity of 
r' = R', 



C^^e-, .-Ce-, (4.11) 

where C is constant and the cylindrical coordinates where used. The unstable direction 
corresponds exactly to the unstable direction in the classical centrifugal instability, given 

by 

I = f 217 (r') + r'f-\ = (/) = ^-{l + Ro), (4.12) 

where r2('r') ~ l/r'" and the expression (|3.17p were used. 

In the magnetic case the Lorentz force enters the balance, which depends on the wave- 
length of the perturbations through the electric currents. Therefore the flows which are 
stable in the non-magnetic case may become unstable in the presence of magnetic field if 
the wavelength of the perturbation can be sufficiently adjusted. The Navier-Stokes and 
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the induction equations for the horizontal rotational MEI take the form 

'di 



(v • V) uo - 20 X V + ihh, 



^ = (b-V)uo 



i/iv, 



(4.13) 



(4.14) 



For e = 1 they correspond to the MRI equations in Balbus & Hawley (1991) with the 
radial wave vector (and thus the vertical components of the perturbations) set to zero (in 
their notation k^, = 0). Of course such restriction makes the considerations of Balbus & 
Hawley (1991) less general. However, it leads to the same stability criterion dVL/dR > 0. 

As previously we can separate the total force on the right hand side of ()4.13p into a 
component along the perturbation vector v, 



fII - 

total 

and a component perpendicular to v 

ih 



ih 



total 



(Vxby 



Vybx) 



2Ro"^w. 



(4.15) 



(4.16) 



Next we find a general form of a j-ih. eigenvector of the matrix in (|2.5p associated with 
the eigenvalue (jj {j = 1,2,3,4), 



vi = -2Ro-^h(jj, 
&J =/i[(£ + l)(e_l-2Ro-i) 



bi ^ -I 



erf {\ + £) + h^ {e+\ + 2Ro"i) - (e + 1) (e^ - (l + 2Ro"i 



2Ro" 



(4.17a) 
(4.176) 
(4.17c) 

(4.17rf) 



Since for any of the four eigenmodcs -b^ — Pc{<^j) ~ 0, where Pc{-) is the characteristic 
polynomial of the matrix in (|2.5p we get (for a single eigenmode) 



total 



-4£Ro"^ [(e + 1) (e - 1 - 2Ro" 



(4.18) 
ai (cf. (l2J|) ) 



We now concentrate on the unstable mode with index 1. Inserting (jj 
into the above expression one can easily show that F^^^^i > if tri > 0, i.e. for > "0+- 
The calculation of the force component perpendicular to the perturbation velocity 
vector for the unstable mode 1 is slightly more involved, but it is straightforward to 
show that 

FLiai = 0- (4.19) 
Thus the mechanism of the instability is essentially the same as in the non-magnetic 
case, i.e. for the instability to be excited there must exist a possibility for a balance of 
all the rotational components of the forces acting on the perturbation, namely strain, 
the Coriolis and Lorcntz forces, so that the perturbation for which such balance occurs 
can be amplified constantly in time by the strain. If such a balance cannot occur, the 
perturbations are amplified and damped in finite periods of time, thus they oscillate and 
do not grow. The unstable direction can be obtained from (j4.17b .b). It can be seen, that in 
the magnetic case the unstable direction depends on the wavelength of the perturbation. 
This is no surprise, since it was already mentioned that since the Lorentz force depends on 
the wavelength through the currents, the balance of the rotational components depends 
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on the choice of k. This explains why an ehiptic vortex, which is stable in the non- 
magnetic case, may become unstable when the magnetic field is turned on with respect 
to perturbations for which the wavelengths must satisfy a certain condition. For the case 
of rotating shear flow, i.e. e = 1, a system with any negative Rossby number (Ro < 0) 
is unstable with respect to horizontal perturbations with the wave number k satisfying 
(cf. (|2.13p ') k^ < — 4RoO^/w^. In a vertically bounded domain, such as an accretion disk, 
there is a lower bound on the allowed values of k, thus such conditions are, in principle, 
modified (cf. Bajer & Mizerski 2011). 

5. Conclusions 

We have studied the properties and physical meaning of the horizontal MERI, i.e. the 
instability of the Euler flow with elliptical streamlines in the presence of external magnetic 
field and background rotation with respect to horizontal perturbations propagating along 
the field lines. The presence of anticyclonic background rotation is necessary for excitation 
of the horizontal instability and we have demonstrated that adding the magnetic field 
changes the range of values of the Rossby number at which the instability develops from 
Ro < —1 in the non-magnetic case to — 47^/(fc^i;^) < Ro < when the magnetic field is 
switched on. 

The main purpose of this paper is the unification of the two, previously unrelated 
instability types. By the use of the local shearing box approximation the direct cor- 
respondence in the stability characteristics was shown between the well known short- 
wavelength MRI and the horizontal rotational MEI in the limit of rotating shear flow. 
The physical conditions under which the direct link between the two instabilities exists 
and the correspondence of the parameters between both instability types were explained 
in detail. The exponent s defining the radial dependence of the angular velocity in the 
disc, ~ l/r*, is related to the local Rossby number in the box by s = — 2Ro and 
the MRI requirement for instability, namely s > corresponds to the rotational MEI 
requirement, i.e. Ro < 0. Moreover, the dispersion relation for rotating shear flow with 
axial magnetic field for the horizontal modes was shown to be exactly the same as for 
the short-wavelength MRI and the well-known feature of the MRI in Keplerian discs, 
that the instability occurs for < kvA/^ < with the maximal growth rate a = 30/4 
achieved for /cua/O = •\/T5/4 is also observed in rotating shear flow. The non-magnetic 
case of the horizontal instability was shown to have the same stability characteristics in 
the limit of rotating shear flow as the local centrifugal instability, with the Rayleigh's 
criterion corresponding to the non-magnetic horizontal instability requirement Ro < — 1. 
Finally, we have also provided a physical interpretation of the horizontal instability, 
based on the balance between the strain generated by the basic elliptic flow and the 
Coriolis force. The unstable and stable directions of the perturbation velocity field as 
functions of Ro and £ were found and the correspondence to the unstable direction of the 
perturbations for localised centrifugal instability was shown. It was also demonstrated 
that instability can only occur, if there is a possibility for the Coriolis force to be balanced 
by a component of the strain to eliminate the rotational tendency in the perturbation 
velocity, which, if present for all times, leads to oscillations. Such balance defines the 
unstable/stable direction and can occur only if > (l + 2Ro~^) , which agrees with 
the results of the linear analysis. 

The full magnetic problem, for which the instability mechanism is essentially the same, 
was also discussed, highlighting the effect of the Lorentz force, which significantly af- 
fects the force balance, making the unstable direction dependent on the wavelength 
of the perturbation. The fact that the Lorentz force can be adjusted by the choice 



14 



K. A. Mizerski and W. Lyra 



of the perturbation wavelength results in a possibility of achieving the balance of the 
rotational components of all three forces, which leads to continuous amplification of 
the perturbation, for negative Rossby numbers even outside the non-magnetic interval 
2/(e — 1) < Ro < — 2/(e + 1). In particular in the rotating shear flow limit, e = 1 all 
systems with Ro < are unstable, hence the presence of magnetic field leads to desta- 
bilisation of all elliptical vortices with Ro < — 1. 

WL gratefully acknowledges partial financial support of the US- American National 
Science Foundation (NSF) under grant no. ASTlO-09802. 
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